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Abstract 

We study a class of subdivision invariant lattice models based on 
the gauge group Z p , with particular emphasis on the four dimensional 
example. This model is based upon the assignment of field variables 
to both the 1- and 2-dimensional simplices of the simplicial complex. 
The property of subdivision invariance is achieved when the coupling 
parameter is quantized and the field configurations are restricted to 
satisfy a type of mod-p flatness condition. By explicit computation of 
the partition function for the manifold i?P 3 x S , we establish that the 
theory has a quantum Hilbert space which differs from the classical 
one. 



ITFA-94-13 
May 1994 - Revised 9/94 



1 Supported by Stichting voor Fundamenteel Onderzoek der Materie (FOM) 
Email: Dannyb@phys.uva.nl 

2 Email: Rakowski@maths.tcd.ie 



1 Introduction 



A series of Z v lattice models was introduced in |T| which had the very spe- 
cial property of being subdivision invariant. This means that the partition 
function is insensitive to successively finer triangulations of the underlying 
simplicial complex. One should regard this property as the discrete analog 
of a continuum quantum field theory being metric independent. The formu- 
lation of these models involved the assignment of field variables to simplices 
of various dimensions. In three dimensions, only link based gauge fields are 
possible and that theory reduced to the abelian Dijkgraaf-Witten model [0]. 
A new four dimensional model was also introduced which involved fields as- 
sociated to both 1- and 2-dimensional simplices of the simplicial complex. 

A crucial element in securing the property of subdivision invariance was 
to restrict the allowed field configurations to those satisfying a certain "flat- 
ness" condition; in addition, a quantization of the coupling parameter was 
also necessary. Solutions to the flatness conditions correspond to simplicial 
cohomology classes of the underlying complex K. The partition function is a 
sum over these classes of a Boltzmann weight which captures a certain kind 
of "intersection" of these field configurations. In our four dimensional model, 
this intersection is between H l (K, Z p ) and H 2 (K, Z p ). 

Our aim here is to develop further the properties of this theory, and 
specifically to establish "non-triviality" in four dimensions. By this we mean 
that we have a topological field theory whose quantum Hilbert space differs 
from the classical one; it is simply a statement about the dependence of 
the theory on the coupling parameter in the Boltzmann weight. This can 
be contrasted with the Dijkgraaf-Witten model with gauge group Z p , where 
there is no distinction between the classical and quantum Hilbert spaces. 
Recall that a topological field theory in d + 1 dimensions associates a Hilbert 
space to each closed ci-manifold. The d + 1 dimensional theory then governs 
the topology changing amplitudes between rf-manifolds which appear on the 
boundary. In ||, such a model was constructed in three dimensions and there 
the dimensions of the quantum Hilbert spaces for various bounding Riemann 
surfaces were related to conformal field theory. The novelty in our models is 
that one can study examples in four and higher dimensions as well. These 
models should also prove useful in the general classification programme of 
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topological field theory. 

After reviewing some general properties, we consider in detail the evalua- 
tion of the partition function on the manifold RP 3 x S 1 which computes the 
dimension of the Hilbert space associated to RP 3 . This can then be compared 
with the simple example of S 3 x S 1 . We relate our pedestrian formulation of 
these theories with the Bockstein operator, and finally we present some of the 
properties associated to 4-manifolds with boundary, including the behaviour 
under connected sum. 

2 General Formalism 

A lattice model is based on a simplicial complex which combinatorially en- 
codes the topological structure of some manifold. Let us recall some of the 
essential ingredients that are required in such a formulation; we refer the 
reader to || |[ U for a more complete account. 

Let V = {vi} denote a finite set of N points which we will refer to as the 
vertices of a simplicial complex. An ordered fc-simplex is an array of k + 1 
distinct vertices which we denote by, 



It will usually be convenient to use simply the indices themselves to label a 
given vertex when no confusion will arise, so the above simplex is denoted 
more economically by [0, - • • , k]. Pictorially, a fc-simplex should be regarded 
as a point, line segment, triangle, or tetrahedron for k equals zero through 
three respectively. A simplex which is spanned by any subset of the vertices 
is called a face of the original simplex. An orientation of a simplex is a choice 
of ordering of its vertices, where we identify orderings that differ by an even 
permutation, but for the models described here we will require an ordering 
of all vertices. One then checks that the invariant we compute is actually 
independent of the choice made in vertex ordering. 

The boundary operator d on the ordered simplex a = [v , ■ ■ ■ ,Vk] is de- 
fined by, 



[vo, ■■■,v k ) 



(1) 



k 




(2) 



i=0 
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where the 'hat' indicates a vertex which has been omitted. It is easy to show 
that the composition of boundary operators is zero; d 2 = 0. 

We model a closed n-dimensional manifold as a collection K = {a,i} of 
n-simplices constructed from the set of vertices V, subject to a few technical 
conditions. Most importantly, every (n — l)-face of any given ra-simplex ap- 
pears as an (n — l)-face of precisely two different n-simplices in the collection 
K. One thinks of the n-simplices then as glued together along (n — l)-faces. 
In order to ensure that the simplicial complex represents a manifold, we re- 
quire the "link" of each vertex to be a combinatorial (n — l)-dimensional 
sphere. We refer the reader to [[| || for a more complete discussion of this 
condition. 

The dynamical variables in the theories we construct will be objects which 
assign an element in the cyclic group Z p = Z/pZ, which we represent as the 
set of integers, 

{0,---,P-l} , (3) 

to ordered simplices of some specified dimension. We call these dynamical 
variables A;-colours with coefficients in Z p , and denote the evaluation of some 
fc-colour B^ on the ordered A;-simplex [0, • • • , k] by 

<B^,[0,---,k]>=B ... k eZ p . (4) 

The superscript (k) will usually be omitted when its value is clear from 
context. It is important to note that we are assigning a Z p element in a way 
which depends on the ordering of vertices in the simplex; we do not have the 
rule = — -B^o , for example. Instead, we shall assume that, 

B§ = -B$j> mod p , (5) 

and similarly extend this to a fc-colour for odd permutations of the vertices. 
The case closest to conventional lattice gauge theory is where a 1-colour 
variable is assigned to every 1-simplex in the complex. 

The coboundary operator 5 acts on the dynamical variables as follows. 
Given a (k — l)-colour, an application of the coboundary operator produces 
an integer in Z, when evaluated on an ordered fc-simplex, namely 

<5B^- 1 \[0,---,k}> = < B,d[0,---,k}> 

= B\23...k — B 23...k + B Q i 3 ...k — ■ ■ ■ ■ (6) 
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We must emphasize that the above sum of integers is not taken with modular 
p arithmetic; it is simply an element in Z. In cases where we will need to take 
some combination mod-p, we will put those terms between square brackets, 
so for example, 

[a + b] = a + b mod p . (7) 

There is also a cup product operation on colours which takes a /c-colour 
and a /-colour and gives an integer in Z when evaluated on a (k + l)- 
ordered simplex: 

< B U C, [0, • • • , k + 1} >= B ... k ■ C k ... k+l . (8) 

Note once again that this product is in Z and the value is not taken mod-p. 

Let us now put these ingredients together and define our theories. First, 
we must be given some oriented simplicial complex K which we take to 
represent a manifold of dimension n. One then has some collection of n- 
simplices defined up to orientation. Take the vertex set of this complex and 
give it an ordering. This is done arbitrarily and we will have to show that our 
construction is independent of this choice, see for example @, || [7j. Now we 
can write down an ordered collection of the n-simplices; each of the simplices 
is written in ascending order and a sign in front of that simplex indicates 
whether that ordering is positively or negatively oriented with respect to the 
orientation of the complex K. Let us denote this ordered set of n-simplices 
by K n , 

K n = J2 , (9) 

i 

where the index i runs over the ordered ra-simplices <jj and e« is a sign which 
indicates the orientation. We will assign a Boltzmann weight W[K n ] to K n 
by taking a product of factors, one for every n-simplex, 

W[K n ] = J] WfaY* . (10) 

i 

Each of the individual factors is a nonzero complex number and will be 
some function of the colours. The details of which colours we use and how 
the function is defined will depend on the particular model. Finally, the 
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partition function, which we will require to be a combinatorial invariant, is 
defined to be a quantity which is proportional to the sum of the Boltzmann 
weights over all colourings, 

z = \Girn £ w[k»\ ■ (n) 

I I colours 

Here \G\ is the order of the gauge group and f(N) is a function of the 
number of simplices of various dimensions. This function will be fixed for 
any given theory by scaling considerations. In the four dimensional model to 
be discussed next, f(N) = N\ where Ni is the number of 1-simplices in the 
simplicial complex. In the three dimensional Dijkgraaf-Witten model (based 
on a single 1-colour field), as formulated in |]|, it is equal to the number of 
vertices N . 



3 State Sum Model in Four Dimensions 

Let us now turn our attention to the four dimensional model of interest. This 
model is based upon the assignment of field variables to both the 1- and 2- 
dimensional simplices of the simplicial complex. The Boltzmann weight of 
an ordered 4-simplex [0, 1, 2, 3, 4] is defined by: 

W[[0, 1, 2, 3, 4]] = exp{/5 < B (2) U <L4 {1) , [0, 1, 2, 3, 4] >} 

= exp{(3B i 2 (A 23 + A 34 - A 24 )} , (12) 

where and A^ are 2- and 1-colour fields, respectively. Here, (3 is a 
complex number which is as yet unrestricted; we shall also find it conve- 
nient to use the scale factor s = exp[/3]. The first item on the agenda is to 
demonstrate that the Boltzmann weight defines a theory which is subdivi- 
sion invariant. As we shall see, this requirement will enforce a quantization 
of the coupling parameter, and lead to a restriction on the allowed colour 
configurations. In order to establish the property of subdivision invariance, 
it is sufficient to show that the Boltzmann weight is invariant under a set of 
moves known as the Alexander moves M . Equivalently, for the case of closed 
manifolds, we can establish invariance by examining the behaviour under a 
set of (k, I) moves [[|, which we now recall. 
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The (k, I) Moves: 



In the four dimensional case of interest here, we have five (k, I) moves, 
with k = l,---,5, and k + I — 6. It suffices to consider the first three 
cases; the (4,2) and (5,1) moves are inverse to the (2,4) and (1,5) moves, 
respectively. 

The (1,5) move: 

This is described by adding a new vertex x to the centre of the 4-simplex 
[0, 1,2,3,4], and linking it to the other 5 vertices. The original 4-simplex is 
then replaced by an assembly of five 4-simplices, written symbolically as: 

[0,1,2,3,4] - [x,l,2,3,4]-[x,0,2,3,4] + [x,0,l,3,4] 

- [rr,0,l,2,4] + M,l,2,3] • (13) 

This move is also known as an Alexander move of type 4. Note also that we 
declare the new vertex x to be the first in the total ordering of all vertices. 

The (2,4) move: 

In this case, two 4-simplices which share a common 3-simplex [0,1,2,3] 
are replaced by four 4-simplices sharing a common 1-simplex [x, y] : 

[a, 0,1, 2,3]- fo, 0,1, 2,3]- 

[x, y, 1, 2, 3] - [x, y, 0, 2, 3] + [x, y, 0, 1, 3] - [x, y, 0, 1, 2] . (14) 

Again, we place the new vertices x, y at the beginning of the vertex list. 

The (3, 3) move: 

[y, z, 0,1,2] -[x,z, 0,1,2] + [x,y, 0,1,2] - 
[x,y,z,l,2]-[x,y,z,0,2] + [x,y,z,0,l] . (15) 

We note that the 2-simplex [0,1,2] is common to the left hand side, with 
[x, y, z] being common to the right. 

For the case of the (1, 5) move, one finds that the Boltzmann weights 
before and after subdivision are related by: 

W[[0, 1, 2, 3, 4]}s- <5Bu5A > M>i,2,3,4]> = W [fa 1; 2 , 3, 4]] (16) 
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W[[x, 0, 2, 3, 4]]" 1 W[[x, 0, 1, 3, 4]] W[[x, 0, 1, 2, 4]]" 1 W[[x, 0, 1, 2, 3]] . 

It is immediately evident that the Boltzmann weight is not generally invariant 
under this move, due to the presence of the added "insertion" on the left 
hand side of ([16]). Our task is therefore to trivialize this unwanted insertion 
factor, and this can indeed be achieved by imposing a restriction on the 
sum over colourings and on the parameter (3. Subdivision invariance of this 
four dimensional theory is now guaranteed by imposing quantization of the 
coupling s p = 1, as well as a restriction of the colourings to those satisfying 
the conditions 

[5B {2) ] = [5A m ] = . (17) 

We shall refer to these restrictions as "flatness" conditions. For example, on 
the 2-simplex [0, 1, 2], we have the restriction on the 1-colour field 

[5A] 012 = [A 12 -A 02 + A 01 ]=0 . (18) 

As a reminder, we note that this particular equation can also be written as 

[A 01 + A l2 ] = A 02 . (19) 

On the 3-simplex [0, 1, 2, 3], the restriction on the 2-colour takes the form: 

[<5£?]oi23 = [-B123 — -B023 + -B013 — -B012] = . (20) 

With these restrictions, the product 5B U 5 A is clearly a multiple of p 2 and 
the above insertion becomes unity. The resulting identity involving the six 
Boltzmann weight factors shall be referred to as the 6W identity. It is worth 
pointing out that invariance is achieved here without the necessity of sum- 
ming over the additional configurations attached to the vertex x. 

It requires little extra work to complete the demonstration of subdivision 
invariance. One first notes that the remaining (k, I) moves also involve six 
Boltzmann weight factors, and it is easy to see that the 6W identity is also 
a statement of invariance under the (2, 4) and (3, 3) moves. 

The subdivision invariant Boltzmann weight for the 4-simplex [0, 1, 2, 3, 4] 
is given by: 

W[[0,l,2,3,4]] =exp{^-B 012 (A 33 + A u - [A 23 + A u ])} , (21) 
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with k G {0,1, 1}. 

At this point, we can reveal that each of the colour fields enjoys a local 
gauge invariance. The gauge transformation of the A field defined on the 
ordered 1-simplex [0, 1] is defined by: 

A' 01 = [A- 5u} 01 = [A 01 + uo} , (22) 

where a; is a 0-colour field defined on the vertices of the complex. For the 
2-colour field B defined on the ordered 2-simplex [0, 1,2], we have a gauge 
transformation given by: 

B'm2 = [B — #A]oi2 = [-B012 — A12 + A 2 — A i] , (23) 

where A is a 1-colour defined on 1-simplices. Our task now is to show that 
the Boltzmann weight for the case of a closed simplicial complex is invariant 
with respect to independent gauge transformations of the A and B fields. 
As we shall see, invariance of the theory under the above transformations is 
not manifest, but requires both the quantization of the coupling parameter, 
together with the restriction on the allowed field configurations. 

Under the transformation of B, one finds that 

S B'USA _ S BUSA S ~5XU8A _ g BUSA s -8(Xu8A) ^4) 

where the first equality uses the fact that 8 A is an integer multiple of p due to 
the flatness constraint, and that s is ap 2 -root of unity. Hence, the Boltzmann 
weight is invariant up to a total boundary term and the product of all these 
cancels for a closed oriented complex. To demonstrate invariance under the 
A field transformation, one first notes the simple identity 

S BU5A = S -8BUA S 8{BUA) _ ^ 

Invariance then follows immediately by the above argument. 

As discussed in the previous section, the Boltzmann weight is initially 
defined for a specific ordering of the vertex set. We recall here a simple 
argument presented in |L(| which can be used to verify that the value of the 
partition function is independent of this choice. 

Let V = {v , ...,t>7v -i} be the vertex set of the complex, I the index set 
I = {0, ...,Nq — 1}, and define a vertex ordering to be a map / : V — > I. 
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Clearly, if /' is a different vertex ordering, then the composition /' o f^ 1 is 
a permutation on the set I. Furthermore, to each permutation there is a 
corresponding vertex ordering. Since any permutation of / can be decom- 
posed as a product of transpositions of consecutive numbers, it suffices to 
show that the Boltzmann weight is invariant when two consecutive values of 
the ordering / are permuted. Our task is therefore to show that the Boltz- 
mann weights defined with an ordering /, and /' = 7r o /, coincide. Here, 
the permutation 7r is defined by — j + 1, ir(j + 1) = j for some j, and 
7r(i) =i if i ^ {j,j + 1}. 

If j and j + 1 label vertices which do not bound a 1-simplex, then the 
Boltzmann weight is clearly invariant. This follows because j and j + 1 are 
simply dummy variables which can be freely exchanged, without affecting 
the orientation of any individual 4-simplex in the complex. 

In order to establish invariance when the vertices labelled j and j + 1 
bound a 1-simplex, we recall the definition of an Alexander move of type 
1. Given an ordered 4-simplex [v ,Vi,V2,v 3 ,V4\, we introduce an additional 
vertex x at the centre of the 1-simplex [i> ,i>i], giving rise to the move 

[v , V!, V 2 , V 3 , V 4 ] -> [X, V!,V 2 , v 3 , v 4 ] - [x, v , v 2 , v 3 , v A ] . (26) 

Since we have shown that the Boltzmann weight is invariant under the (k, I) 
moves, it is equivalently invariant under all Alexander moves. Thus we are 
free to perform an Alexander move of type 1 on the 1-simplex with vertices 
labelled by j and j + 1. This has the effect that these vertices no longer 
bound a 1-simplex, and by the above argument j and j + 1 can then be 
interchanged leaving the Boltzmann weight invariant. In order to recover the 
original complex with the permuted vertex ordering, one simply performs the 
inverse Alexander move of type 1. 

We have already shown that the Boltzmann weight is invariant under all 
the (k, I) subdivision moves. However, recall that to achieve subdivision in- 
variance, we are required to restrict the allowed field configurations to those 
satisfying the appropriate flatness conditions. This is effected in the state 
sum through the insertion of a set of delta functions which implement the 
required restrictions. It remains to check the behaviour of these delta func- 
tions under the (k, I) moves. As we shall see, the true subdivision invariant 
partition function is given by including a certain scaling factor, as discussed 



9 



in relation to equation (|Tl|). This takes into account the redundancy in the 
assembly of delta functions which are present under subdivision. 

In order to determine the correct scaling factor, we need to examine the 
behaviour of both the A and B delta functions with respect to the (k, I) 
moves. If we denote by AiVj the increase in the number of z-simplices due to 
a (k, I) move, then it is straightforward to check that under the (1, 5) move 
we have: 

AiVo = 1 
AiVx = 5 
AiV 2 = 10 
AiV 3 = 10 

AiV 4 = 4 . (27) 
The changes under the (2, 4) move are given by 

AiV = 
AiVi = 1 
AiV 2 = 4 
A7V 3 = 5 

AN A = 2 , (28) 

and of course under the (3, 3) move we have AiVj = 0, for all i. 

Let us now consider the behaviour of the B delta functions under subdi- 
vision. We will first collect some formulas and then put the results together 
to determine the form of the scaling factor f(N) referred to in equation (ffT|). 
If we denote the additional ten B fields present after a (1, 5) move by: 

I — {B x oi, B x02 , B x03 , B x0 4, B xl2 , B xl3 , B x i4, B x2 3, B x2 ^ B x34 } , (29) 

then one readily finds that summation over these fields yields the result 

1 ^E 5 ([ 5S ]-012)5([5S] :r oi3)5([55] a;014 )5([5S] x02 3)5([5 J B] :c02 4) 

M i 

*([^]^) < j([ < js] al2 3)«y([«ys] a!l2 4)(j([«y5] ail 34)(y([«5 J B] a!a 34) 

5([5B] 0123 )5([5B] 0124 )5([5B] 0rM )5([5B] 02U )5([5B} 12U ) 
= 5([5B] 0123 )5([5B] ou ^B]oi 3 M{5B] 02U )5([5B] 12U ) . (30) 
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Here, the assembly of delta functions on the right and left hand sides above 
represent the situation before and after subdivision. We specify that the 
modulo-p delta function is defined by 

8(x) = { 1 if X = ° mod ' p (31) 
^ ' [ otherwise. ^ ' 

For the case of the (2,4) move, one finds that summation over the addi- 
tional four B fields 

I — {B xy o, B xy i, B xy 2, B xy3 } , (32) 

produces the result: 

^^^[^^oi)^^]^)^^]^)^^]^)^^]^) 
*([^]^)(J([(JS] a! oi2)*([^] a; o 1 3)<J([(JS] a! o23)5([5S] !rl2 3) 

5([5B]y 012 )5([8B} y(m )8([8B}y023)5([8B} yl23) 

^8([8B} x01 2)8([8B} x013 )8([8B} x0 2 3 )8([8B] xl 2 3 )8([8B}y 01 2) 

5([5B] y013 )5([5B] y023 )5([5B] yl 2 3 )5([5B] 0123 ) . (33) 

Turning now to the delta function insertions for the A field, we proceed 
in a similar manner. The additional A fields present after a (1, 5) move are: 

/ = {A x0 , A xl , A x2 , A x3 , A x4 } . (34) 

One verifies that the following relation holds: 

^E*(M]xOl)«([M]*(B)«(M]*Os)«5(M]«M)*(M]*12) 
H I 

<J(M] iel 3)<J([<JA] :Bl4 )<J([^] ie2 3)(J([(JA] !B2 4)(J([(JA] x 34) 
8([5A} 012 )5([5A} 013 )5([5A}ou)8([5A} 023 )5([5A}o24) 

5([5A] 03 M[5A] 123 )5([5A] 124 )5([^]nM[5A]2u) 

= (J(M]oi 2 )<J([<JA]oi 3 )<J(M]oi4)(J([^]o23)5([^]o24) 

(J(M]o34)5(M] 123 )(J([(JA] 124 )(J([(JA] 1 34)(J(M] 2 34) • (35) 
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Finally, we treat the (2,4) move for the A field. There is a single addi- 
tional A field I = {A xy } which is present after subdivision. Summation over 
this field produces the result: 

^5(M]^ )5(M] XJ/1 )^M]^ 2 )5(M]^3)5(M]012) 
5([M]0l3)5([M]023)5(M]l23)5(M],0l)5(M,02) 

5([5A] x03 )5([5A} xl2 )5([5A} xl3 )5([5A} x23 )5([5A} y01 ) 
5([5A] y02 )5([5A} y03 )S([SA] yl2 )5([SA] yl3 )5([5A} y23 ) 

= 5([M] 012 )5([M]013)5(M]023)5(M]l23)5(M,0l) 

5([SA] x02 )S([SA] x03 )S([SA] xl2 )5([5A] xl3 )5([5A} x23 ) 
6([6A} y01 )6([6A} y02 )6([6A} y03 )6([6A} yU )5([5A} yl3 ) 
S([SA} y23 ) . (36) 

We can now establish the correctly scaled subdivision invariant partition 
function by combining the previous results. Under the (1,5) move, we see 
that a factor |G| 5 must be accounted for in the combined A and B sectors, 
and a factor of |G| X under the (2,4) move. But this is precisely how the 
number of 1-simplices changes under these moves. If the partition function 
of (|Tl|) is chosen to have f(N) = N±, then it defines a subdivision invariant 
quantity. Specifically, we have 

W[K n ) , (37) 

M flat 

where we denote the set of allowed colours satisfying the flatness conditions 
by flat. Clearly, at the trivial s — 1 root of unity (k = in equation (PTj)), 
the value of the partition function simply counts the number of solutions to 
the flatness conditions. Our main goal is in achieving interesting behaviour 
at the non-trivial roots of unity where different phase factors can occur. 

4 Evaluation of the Partition Function 

The models described in the preceding sections require that a space be pre- 
sented as a simplicial complex for their formulation. It is clear that one 
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can only hope for a non-trivial partition function - one in which the phases 
(Boltzmann weights) are not all unity - when the field configurations we sum 
over are sufficiently interesting. This means that we need solutions to the 
flatness equations ( |17D which would not be solutions in the "strong sense" if 
the mod-p brackets had been removed. Perhaps the simplest example in four 
dimensions is the space RP 3 x S 1 , and we will give here a rather detailed 
exposition of its simplicial description. 

Let us begin by presenting an economical simplicial complex for the man- 
ifold RP 3 . A complex with a minimal number of 11 vertices has been given 
in and we label its vertices by elements in the set {0, 1, 9, a}. The 
complex is fully determined by specifying the 3-simplices; these are 40 in 
number and are given explicitly by, 



+ [0,2, 9, a] + [0,2,3,9] - 


[0. 


2. 


3,7] 




[0, 2, 7, a] + 


[0, 5, 7, a] 


-[0,4, 5, 7] + [0,1, 4, 5] + 


[0. 


1, 


3,4] 




[0,1,3,9] + 


[0,1,6,9] 


+ [0,1,5,6] - [0,5, 6, a] - 


[o 


6 


9, a 


+ 


[4, 6, 9, a] + 


[4,6,7,9] 


-[4, 5, 7, 9] + [5, 7,8,9] - 


[5. 


7, 


8, a] 


+ 


[1, 7, 8, al- 


[1,7,8,9] 


-[1,6, 7, 9] + [1,2, 6, 7] + 


[1- 


2, 


5,6] 


+ 


ii, 2,4,5] - 


[1,2, 4, a] 


+ [1,3, 4, a] - [1,3, 8, a] + 


[1 


3, 


8,9; 


+ 


[3,5,8,9]- 


2 . 3 j 5 ; 9 


+ [2,4,5,9] + [2, 4, 9, a] - 


[3. 


5, 


8, a] 


+ 


[3, 5, 6, a] — 


[3, 4, 6, a] 


+ [3,4,6,7] + [2,3,6,7] + 


[2. 


3. 


5,6] 


+ 


[1,2, 7, a]- 


[0,3,4, 7] 



(38) 



where the signs denote the relative orientations of each simplex. Of course, 
the lower dimensional simplices are given by all those which appear as sub- 
simplices in the above list. This complex contains 51 1-simplices and 80 
2-simplices in addition to the 11 vertices and 40 3-simplices already tabu- 
lated. The Euler number is zero as required for a closed 3-manifold. One 
also easily checks that the boundary of the above complex vanishes and that 
these 3-simplices are glued together along paired 2-simplices. 

Constructing the complex for RP 3 x S 1 is straightforward. We begin by 
imagining the above complex of 40 3-simplices displayed horizontally. To 
each of those we add a new vertex and construct a vertical tower beneath 
which contains a total stack of 12 4-simplices; this is the S 1 direction which 
gets glued to the top along a common 3-simplex. So, for example, the tower 
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beneath +[0, 1, 3, 4] is given explicitly by, 



+ [0 . 


1 , 


3 , 


4 , 


0'] 


-[1 


3 . 


4 • 


0' 


i'; 


+ [3 . 


4 ,0' 


1' 


3' 


"[4 


0' 


l' 


3' 


4' 


+[0' 


1' 


3' 


4' 


,0" 


-[!' 


3' 


4' 


,0" 


1" 


+[3' 


4' 


0" 


1" 


3"] 




0" 


1" 


3" 


4" 


+ [0". 


1", 


3", 


4", 


0] 


-[1" 


3". 


4". 


, 


1 ] 


+[3", 


4", 


, 


1 , 


3] 


-[4" 


. 


1 • 


3 , 


4] 



(39) 



One sees that for each vertex x in the original complex for RP 3 , two new 
vertices x' and x" are required in this presentation of RP 3 x S 1 . The total 
vertex set now is then, 



{0, ...,a,0', ...,a',0", ...,a"} 



(40) 



and contains 33 elements. It is straightforward, though tedious, to enumer- 
ate all simplices in this complex. The number of each simplex type in this 
simplicial complex for RP 3 x S* 1 is, 

— simplices 33 (41) 

1 — simplices 339 

2 — simplices 1026 

3 — simplices 1200 

4 — simplices 480 . 

In order to compute the partition function in these theories, we need to 
first determine the admissible field configurations. This means finding the 
gauge inequivalent solutions to the equations, 



[5 A] 012 = [A 12 - A 02 + A i] = 

[8 B\ 0123 — [B\23 — Bq23 + A)13 — -B012] 



(42) 
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where gauge transformations are given by, 

A' 01 = [A 01 - out + uj ] (43) 

-^012 = [-^012 — ^12 + ^02 — Aoi] • 

We remind the reader that the brackets in these equations denote that the 
quantity inside is to be taken mod-p. 

The number of gauge inequivalent solutions to these equations is in cor- 
respondence with the first and second cohomology groups of the complex K 
with coefficients in Z p . Beginning with the well known homology groups with 
integer coefficients for RP 3 and S* 1 [||, f| , 

H (RP 3 ) = H 3 (RP 3 ) = HoiS 1 ) = H^S 1 ) = Z , 

H^RP 3 ) = Z 2 , H 2 (RP 3 ) = , (44) 

the Eilenberg-Zilber theorem, 

H n (XxY)= HiWQHAY)® Yl Tor(H p (X),H q (Y)) , (45) 

i+j=n p+q=n—l 

computes the homology of the product K = RP 3 x S 1 , and one finds 

H Q (K) = H 3 (K) = H 4 (K) = Z , H 1 (K) = Z®Z 2 , H 2 (K) = Z 2 . (46) 

The universal coefficient theorem for cohomology, 

H n {X, G) = Hom(H n (X), G) © Ext{H n ^{X), G) , (47) 

then gives the required cohomology groups which for G = Z p are, 

H°(K,Z p )=H\K,Z p )= Z p (48) 

H\K, Z p ) = H 3 (K, Z p ) = I l P r, 1° IP ° dd 
v ' pJ K ' v> \Z P ®Z 2 for p even 



H 2 (K,Z p ) 



for p odd 

Z 2 © Z 2 for p even. 



Now, the partition function in this four dimensional model essentially 
amounts to a sum over field configurations which represent inequivalent 
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classes in the cohomology groups H 1 ^, Z p ) and H 2 (K, Z p ). While the scale 
factor we have introduced, l/G Nl , defines a subdivision invariant quantity, 
one could also adopt a different normalization where the partition function is 
precisely proportional - in a way independent of the simplicial complex - to a 
sum over these classes. To relate our original partition function to the latter, 
we need to count carefully gauge equivalent copies of all field configurations. 
Also, knowing the number of allowed gauge constraints is important for the 
purposes of explicitly finding all solutions. 

For link based fields, the counting of gauge copies is the same as in lattice 
gauge theory where a different copy of the gauge group is assigned to each 
vertex. The gauge transformation here is, 

A' = [A-5u} , (49) 

and what we seek is the dimension of the image of the map 

5° : C°(K) -> C\K) . (50) 

Here we have explicitly attached a superscript to 5 to denote the restriction 
to C°. But 5°(C°) is isomorphic to, 

C°/Ker(8°) (51) 

and the kernel of 5° is 1 dimensional for a connected complex (same as 
H°(K, Z p )). One then sees that the image of this map has dimension N — l, 
and this is the number of links we can gauge fix. 

For the 2-simplex field B, the counting is only slightly more difficult. As 
a gauge field, we are assigning an element in the gauge group to each link, 
and we seek the dimension of the image of the map, 

5 1 : C\K) -> C 2 (K) , (52) 

which will tell us how many of the 2-simplex fields can be gauged away. Let 
us restrict the following discussion to the case of p a prime number so Z p is a 
field. The kernel of 5 1 , the 1-cocycles, is then parametrized by the image of 
5° together with H l (K, Z p ). When p is prime, the later cohomology group 
is then a sum of copies of Z p ; let h 1 denote the number of these copies, i.e. 
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the dimension of H 1 (K,Z P ) as a vector space over Z p . Hence the image of 
S 1 has dimension, 

Ni - (N - 1) - h 1 . (53) 

Putting these numbers for the maximal trees together, we can, for |G| = p a 
prime number, write the partition function ([37]) as 

Z = jA^ E WW , (54) 

M flat 1 

where flat' indicates the sum is over the gauge inequivalent configurations 
which are the cohomology classes. 

In the case of the complex for RP 3 x S 1 , we have 1026 equations for 
the A field and 1200 for B. These are highly redundant due to the Bianchi 
identities, but nevertheless, the number is quite large and one needs to make 
maximal use of the gauge freedom. For the link based gauge field A, one is 
allowed to set to zero (gauge fix) the fields on a maximal tree. A maximal 
tree is any maximal set of links which contains no closed loops, and as we 
saw above, that number is always one less than the number of vertices. It is 
trivial to pick such a set by inspection, and in this case we can gauge fix the 
A field on 32 1-simplices. 

For the 2-colour field B, the situation is somewhat more intricate. In 
practice, it is not easy to identify such a set by inspection of the complex. 
Instead, we solved this problem by associating a vector of length 339 (one 
place for each link field A) to each of the 1026 2-simplices; this vector then 
represents a gauge transformation. By using a program in Mathematica 
[ |i~2]j , we could find a maximal number of linearly independent vectors which 



was found to be 306 for our complex. There is one additional complication 
however. The RowReduce routine in Mathematica gives vectors which are 
linearly independent over the real numbers, and we seek a set which is linearly 
independent over Z p . One indeed finds a single vector in that set which is 
not linearly independent for all p. Our gauge choice then amounts to 305 
conditions. This is the number given by (|53"D for our complex ([0]) when 
p = 2. While it is not a maximal tree in general, it is an allowed choice for 
all p. 

Another check on the gauge choice we have made here is to nominate 
one link variable as the independent one for each of the gauge conditions we 
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seek to impose. One then shows that the 305 choices found previously can 
be made with no duplications. 

At this stage, solving the equations subject to a maximal gauge choice 
is not difficult, though it is somewhat tedious. Typically, repeated use of 
the gauge conditions forced most other fields to vanish. For the A field, we 
found that the nonzero pieces could be parameterized in terms of two mod-p 
variables a and x, where [2 a] = and x was unconstrained. The same is 
also true for the B field, where we parametrize the solution in terms of b and 
y, with [2 b] = and y unconstrained. 

For each of these field configurations, one then computes the Boltzmann 
weight which is a product of 480 factors, one for each 4-simplex in the com- 
plex. We find the Boltzmann weight 

,2irik n , , . . 

exp [—— 2 6a . (55) 

One notices immediately that the Boltzmann weight is independent of the 
x and y parameters in the general solution. Since the p odd case has no 
non-zero solutions for a and b, there are no non-trivial phases. For p even, 
the sum over a and b yields 

3 + (-l) k = 2 • 2 5 ^ , (56) 

and the partition function at s = exp[2iiik / p 2 } is given by, 

The symbol S p (k) denotes the mod-p delta function; its value is 1 if k — 
mod-p, and otherwise. In detail the calculation for p even takes the form, 

Z = ^ 30 VV(3 + (-l) fe ) • (58) 

The number 339 comes from the number of 1-simplices in the complex, the 
factors with 305 and 32 take into account the gauge equivalent copies of 
the solutions to the flatness equations, and the remaining factors come from 
summing the solutions over a, x, b, and y. 
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One can compare this result to that obtained for the 4-sphere S 4 . In this 
case, a complex is easily given as the boundary of a single 5-simplex; one has 
the following data, 



— simplices 

1 — simplices 

2 — simplices 

3 — simplices 

4 — simplices 



6 



(59) 
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20 



15 



6 



The calculation of the partition function is easily seen to take the form, 



though there are no interesting solutions to the flatness equations and hence 
no possibility of non-trivial phases in this case. Hence, the value of the 
partition function is independent of k. 

One triangulates the space S 3 x S 1 in the same manner as for the projec- 
tive space we have already considered. There is no possibility of any phases 
and a straightforward calculation gives, 



It is also straightforward to carry out calculations on the spaces L(p, q) x 
S 1 , where L(p,q) is a lens space |J, though the triangulations |ll| get pro- 
gressively larger. In this series, RP 3 appears as L(2, 1). We have also done 
the analogous computation of the partition function for L(5, 1) x S 1 . We find 
that the Boltzmann weight depends on two variables a and 6 which must sat- 
isfy [5a] = [56] = 0, so there are no non-unit phases when p is not a multiple 
of 5. When p is a multiple of 5, the Boltzmann weight takes the form 



z[5 4 ] = iA 5 i = i 



(60) 



Z[S 3 x5'] = l 



(61) 



exp[ 



2irik 



5 6 a] 



(62) 



The partition function becomes, 




. 5^5 (fc) f or p a multiple of 5 
otherwise. 



(63) 
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From these examples, we see that there is generally a dependence of the 
partition function on the coupling parameter. Since the partition function 

on M, x S 1 
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we see that 



/ 3 a j gives the dimension of the Hilbert space 
the quantum Hilbert space associated to M 3 differs from the classical one 
(* = 0). 



5 The Bockstein Operator 

At first sight, the construction of this general class of models that we are 
considering in this paper may seem mathematically unorthodox. The mo- 
tivation stemmed from a desire to realize discrete Chern-Simons and BF 
theories from a concrete point of view. The intuition was that one should 
make use of the coboundary operator on simplicial cochains in some fashion, 
but we did so with a cup product which differed from the usual one in so far 
as we did not take the product mod-p. However, we can, in retrospect, make 
an observation which brings the whole construction into orthodoxy, and this 



is the connection with the Bockstein operator fll5| ; the homotopy-type nature 
of these models is then transparent. We shall restrict attention here to the 
case of closed manifolds, so that the fundamental class exists. 

Let x be an element of the simplicial cohomology group H q (K, Z p ), and 
let x G C q (K, Z) denote a representative of integral cochain. Since 

[5 x] = 0, this means 5x = pu, for some integral (q + l)-cochain it. The 
Bockstein operator 

P:H q (K,Z p )^H q+1 (K,Z p ) (64) 

is defined by 

P(x) = [-6x] = [it] . (65) 
V 

In terms of the Bockstein operator and normal cup product, we can then 
rewrite the Boltzmann weight (EH]) quite simply as, 

exp[ < B U f3(A), a A >] , (66) 

V 
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where 04 is a 4-simplex. Of course, the Boltzmann weight for the abelian 
Dijkgraaf-Witten |2| theory can also be so written; 

exp[^ < A U 13(A), a 3 >] , (67) 

where the connection with Chern-Simons theory is striking. The key observa- 
tion is that the Bockstein operator and not the simple coboundary operator 
is what is relevant in the construction of these models with gauge group Z p . 

The extension of this theory to 5 and higher dimensions is then transpar- 
ent, and one takes the Boltzmann weight, 

exp[— — < AU(3(A)U---U(3(A), a 2m+ i >} , (68) 

where we have m factors of (3(A). It is worth noting that theories whose par- 
tition functions lead to Gauss sums such as the 3d Dijkgraaf-Witten theory 
will also appear in 7d where the link field A in (|67|) becomes a 3-colour field. 
In this case, partition functions will generally be complex valued unlike the 
the class of B U (3(A) theories where (for p prime at least) they are real. 



6 Manifolds with Boundary 

A general axiomatic framework was presented for topological quantum field 
theory (TQFT) in |13|]; see also |b|, [Tj| [T7|]. As we have seen, the general 
class of abelian models considered here can be expressed in terms of standard 
modulo-p cohomological operations, and one expects the axioms of TQFT 
to be satisfied by these models. 

Let K be a 4-manifold with boundary OK. The partition function of 
these models is well defined, and represents a transition amplitude when we 
specify the field configurations on the boundary components. We take the 
field configurations on the boundary to be flat, and define the allowed field 
configuration on K to be all flat configurations which extend those specified 
on the boundary. The partition function remains subdivision invariant as 



long as we keep the triangulation on the boundary fixed ||i7 |. It is convenient, 
however, to rescale the partition function by 

Z'[K) = w NlidK) Z[K) , (69) 
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where w = J\G\, and Ni(dK) is the number of 1-simplices on the boundary 
OK. This scaling gives the following gluing rule, 



Z'[K,r 1 ,r 2 ]=J2 Z\K u t x M-Z\K^M . (70) 

Here K\ is a cobordism between boundary manifolds Si and £ 3 with fixed 
flat field configurations T\ and r 3 , and similarly for K 2 . K represents a 
composition of K\ and K 2 and the above sum is over all intermediate flat 
field configurations. 

Consider now a gauge transformation of the B field, when a boundary 
dK is present. From ([24]), we see that the Boltzmann weights are related by 
a phase factor depending only on the boundary values of the fields, namely: 

S <B'USA,K> _ s <Bu6A,K> s -<\u6A,dK> ^1) 

When computing the partition function on K, we sum over all allowed field 
configurations with fixed boundary data, and thus we see that the parti- 
tion function also transforms with this phase factor. It is equally simple to 
determine the behaviour under a gauge transformation of the A field. 

For the purposes of illustration let us consider the case of a 4-manifold 
of the form M 3 x /, where I is the unit interval, and M 3 is some bounding 
3-manifold. The value of the partition function then gives a transition ampli- 
tude between the two copies of M 3 . Given that H*(M 3 x J, Z p ) = H*(M 3 , Z p ), 
we know that the transition matrix Z-j[M 3 x J] must be diagonal. Moreover, 
because of subdivision invariance, 

Z' lf [M 3 xI] = J2 Z ^ M * x J ] ■ ^/t M 3 x /] , (72) 

3 

which shows that any diagonal element can only be or 1. Since the two 
copies of the bounding manifold appear with opposite orientation, the be- 
haviour of the partition function under changes of cohomology representatives 
(gauge transformations) is given by 

Z' v j, = exp[i«i] Z[ j exp[— iotf] , (73) 

where Z-f denotes the transition amplitude between the initial and final copy 
of M 3 . It is then clear, for example, that the diagonal elements along with 
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the trace of are gauge invariant quantities. In particular, we have the 
result 



which we can interpret as the dimension of the Hilbert space associated to M 3 . 
In taking the above trace, one must take into account the gauge equivalent 
copies of the fixed boundary data. 

Indeed, the value of the partition function for the manifold M 3 x S 1 can 
be obtained more easily by first performing the computation on M 3 x /, and 
then taking the trace. In this way, one sees that a vertical tower (see (p9|)) 
of only four 4-simplices is required. 

Up to gauge equivalence, field configurations on any single boundary com- 
ponent M 3 are in one to one correspondence with the set H 1 (M 3 , Z p ) x 
H 2 (M 3 , Z p ). We can define a complex vector space V(M 3 ) associated to M 3 
by taking it to be the vector space freely generated by this set of field config- 
urations. This we might call the classical Hilbert space of M 3 0. However, 
the partition function on the cylinder represents a transition amplitude, and 
the map Z[^\M 3 x I] may well have a non-zero kernel. The quantum Hilbert 
space H(M 3 ) is defined to be, 



and its dimension is given by the above trace (|74|) . Thus, our computations 
show that the quantum Hilbert space of these models is generically different 
from the classical Hilbert space. By classical Hilbert space one simply refers 
to the situation where k = and the kernel of Z[j on the cylinder always 
vanishes. 

It is interesting to go further and identify precisely the zero modes in the 
examples we have computed. For the p = 2, k = 1 theory with RP 3 bound- 
aries, the calculation of the partition function on RP 3 x S 1 indicates a Hilbert 
space of dimension 2, compared to the classical k = result of 4. Hence there 
must be two zero modes in the "propagator" on the cylinder. One finds that 
these zero modes correspond to the non-trivial H 1 (RP 3 , Z2) configuration 
which means that the Hilbert space is in correspondence with H 2 (RP 3 , Z 2 ). 
Whether this is a general phenomenon, or is something peculiar to the lens 
spaces, is not known. 



£ Z' H [M 3 x /] = Z[M 3 x S 1 ] , 



(74) 



H(M 3 ) = V(M 3 )/Ker(Z'[M 3 x /]) , 



(75) 
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Let us now examine the behaviour of the partition function with respect 
to the connected sum of manifolds. For manifolds M\ and M 2 , the connected 
sum is denoted by M = M\j^M 2 . The manifold M is produced by first 
excising a 4-ball from each of the components Mi and M 2 , which are then 
identified along their common S 3 boundary. According to the tenets of the 
axiomatic approach, the partition function on a closed manifold can be com- 
puted by cutting the manifold along a common boundary, and then taking 
the pairing between the state vectors in the dual Hilbert spaces. 

Now, in order to obtain a relationship between the partition function 
Z[M] and that of its components Z[M\] and Z[M 2 ], we recall the value of 
the partition function Z[S 3 x S 1 ] = 1. We thus see that the S 3 Hilbert space 
is 1-dimensional. It is then a simple consequence of 1-dimensional linear 



algebra |I8] to see that the following relationship holds: 

Z\M X #M 2 \ Z[S 4 ] = Z[Mi] Z[M 2 ] . (76) 

Since Z[S 4 ] = 1 for the four dimensional model under consideration, the 
partition function behaves multiplicatively under connected sum. 



7 Concluding Remarks 



The four dimensional model we have considered here is part of a generic 
construction available in all dimensions HJ. The cornerstone of these state 
sum models is a partition function which is a sum over simplicial cohomology 
classes of a certain Boltzmann weight. This phase factor arises from the 
modulo-p valued intersection form between those classes. In particular, we 
have seen that the relevant "kinetic" operator used to define these actions 
is provided by the Bockstein coboundary operator. We remark that one can 



also consider the "non-kinetic" type models [I5| as providing observables for 
the theories presented here. In three dimensions, it reduces to the Dijkgraaf- 
Witten model || which is related to group cohomology. Such a connection 
is not transparent in general. 

We have shown here that the four dimensional model is indeed non-trivial 
in the sense that interesting phases can be obtained; without them the model 
only counts cohomology classes. In this regard, one sees that the dimensions 
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of the quantum Hilbert spaces are in general different from the classical 
dimensions. By construction, these models lead to piecewise linear invariants; 
however, with the insight that they can be formulated (in the closed case) in 
terms of the Bockstein operator, one sees that the models yield homotopy- 
type invariants. 

In [0, a four dimensional subdivision invariant model was described in 



terms of combinatorial data. Subsequently, it was established in |20], |2T[ , that 
the partition function was expressible in terms of the Euler and Pontryagin 
numbers, and as such encoded classical topological data. On the other hand, 
this model could then be viewed as describing classical invariants in terms of 
a quantum state sum. The models presented here can similarly be viewed as 
providing a quantum state sum formulation of classical modulo-p cohomolog- 
ical data. Perhaps it is also worth remarking on how these models differ from 
the structures presented in Quite apart from having to address issues 
of regularizing the formally divergent path integrals of those models, one is 
also dealing with cohomology with real coefficients; as such, the models are 
insensitive to the presence of any torsion subgroups. However, as we have 
seen for the models discussed here, the essence of non-triviality lies in the 
presence of torsion in the cohomology groups. 
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